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Abstract 

We calculate the form factors for the semileptonic decays of heavy-light pseudoscalar mesons 
in partially quenched staggered chiral perturbation theory (SXPT), working to leading order in 
l/mg, where mq is the heavy quark mass. We take the light meson in the final state to be 
a pseudoscalar corresponding to the exact chiral symmetry of staggered quarks. The treatment 
assumes the validity of the standard prescription for representing the staggered "fourth root trick" 
within SXPT by insertions of factors of 1/4 for each sea quark loop. Our calculation is based on 
an existing partially quenched continuum chiral perturbation theory calculation with degenerate 
sea quarks by Becirevic, Prelovsek and Zupan, which we generalize to the staggered (and non- 
degenerate) case. As a by-product, we obtain the continuum partially quenched results with 
non-degenerate sea quarks. We analyze the effects of non-leading chiral terms, and find a relation 
among the coefficients governing the analytic valence mass dependence at this order. Our results 
are useful in analyzing lattice computations of form factors i? — > vr and D — > AT when the light 
quarks are simulated with the staggered action. 

PACS numbers: 12.39.Fe,12.39.Hg, ll.30.Rd, 12.38.Gc 



I. INTRODUCTION 



Extraction of the CKM matrix elements \ Vub\ and \ Vcs \ from the experimentally measured 
semileptonic decay rates for B nil' and D Klv requires reliable theoretical calculations 
of the corresponding hadronic matrix elements. Recently, there has been significant progress 
in computing these matrix elements on the lattice, with good control of the systematic 
uncertainties Since computation time increases as a high power of the inverse 

quark mass, the light (m, d) quark masses used in the simulations are heavier than in nature, 
and a chiral extrapolation is necessary to obtain physical results. To keep systematic errors 
small, the simulated m, d masses should be well into the chiral regime, giving pion masses 
~ 300 MeV or lighter. Such masses in lattice calculations of leptonic and semileptonic heavy- 
light decays are accessible with staggered quarks I^, [?, JO, Q. The trade-off for this 
benefit is the fact that staggered quarks do not fully remove the species doubling that occurs 
for lattice fermions; for every flavor of lattice quark, there are four "tastes," which are related 
in the continuum by an SU{A) symmetry (or an SU{A)l x SU{A)b. symmetry in the massless 
case). The taste symmetry is broken at non-zero lattice spacing a by terms of order a^. 

The breaking of taste symmetry on the lattice implies that one must take into account 



taste- violations in the chiral extrapolations, leading to a joint extrapolation in bot 
masses and the lattice spacing. Staggered chiral perturbation theory (SXPT) 
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le quark 



allows us to make such extrapolations systematic. For quantities with 

.Q, 16 



leavy quarks, one 



17, bJ into SXPT. 



must also incorporate Heavy Quark Effective Theory (HQET) 
This has been done in Ref. |l9l|. and then applied to leptonic heavy-light decays. Here, we 
extend the analysis of Ref. |l9| to the semileptonic case. 

In addition to the practical implications of taste symmetry violations for chiral extrap- 
olations, the violations lead to a potentially more serious theoretical concern. Simulations 
such as Refs. js], 6, Q, [s], lo], |lo|, Q take the fourth root of the staggered quark determinant 
2l| in an attempt to obtain a single taste per quark flavor in the continuum limit. Were the 
taste symmetry exact at finite lattice spacing, the fourth root prescription would obviously 
accomplish the desired goal, since it would be equivalent to using a local Dirac operator 
obtained by projecting the staggered operator onto a single-taste subspace. Because the 
taste symmetry is broken, however, the fourth root is necessarily a nonlocal operation at 
non-zero lattice spacing 221]. The question of whether the rooted theory is in the correct 
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ass therefore becomes nontrivial. Nevertheless, there are strong theoretical 
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261 ] in the interacting theory, as well as free-theory and numerical 



that the fourth-root trick is valid, i.e., that it produces QCD in the continuum 
The current paper does not actually need to assume that the rooting procedure itself is 



valid. ^ Instead, like previous SXPT calculations for the rooted theory [l2|, [l3|, [l9|, l28|, it 
requires a narrower assumption: that the rooting can be represented at the chiral level by 
multiplying each sea quark loop by a factor of 1/4. This can be accornplished by a quark 
flow analysis 29|, or, more systematically, by use of the replica trick 30|. In Ref. 2J], it 
was shown that the correctness of this representation of the fourth root in SXPT follows in 
turn from certain — in our opinion, rather plausible — assumptions. As such, we assume 
here that this representation is valid. Fitting lattice quantities to SXPT formulae (as in 
Refs. [lo|, Q) provides an additional empirical test of the validity of this representation. 

The main purpose of the current paper is to find SXPT expressions for the form factors of 
the semileptonic decay B Piu, where P is some light pseudoscalar meson, which we will 
refer to generically as a "pion." We consider first the partially quenched case, and obtain the 
full QCD results afterward by taking the limit where valence masses equal the sea masses. 
The 5 is a heavy-light meson made up of a 6 heavy quark and a valence light quark spectator 
of flavor X] we use the notation when confusion as to the identity of the light spectator 
could arise. The P meson (more precisely Pxy) is composed of two light valence quarks, of 
flavor X and y. For simplicity we consider only the case where the outgoing pion is (flavor) 
charged; in other words x ^ y. The flavor structure of the weak current responsible for the 
decay is 1/7^6. 

In our calculation, we take the heavy quark mass mq to be large compared to Aqcd and 
work to leading order in the I/ttlq expansion. Our analysis also applies when the heavy 
quark is a c (z.e., to -D mesons), but we use B to denote the heavy-light meson to stress 
the fact that only lowest order terms HQET are kept. For D mesons, of course, the higher 
order terms omitted here would be more important than for B mesons. 

Discretization errors coming from the heavy quark are not included in the current calcu- 
lations. We assume that such errors will be estimated independently, using HQET as the 



^ Of course, were the fourth root trick to prove invahd, the motivation for the current work would be lost. 
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effective-theory description of the lattice heavy quark 3l|] . It is expected that the errors from 
staggered quark taste- violations, which are considered here, are significantly more important 
at most currently accessible lattice spacings than the heavy-quark errors ^. However, since 
taste-violations decrease rapidly^ when the lattice spacing is reduced, this may change in 
the not too distant future. In any case, the precise quantification of the total discretization 
error will always require simulation at several lattice spacings. 

An additional practical constraint on the current calculation is that antQ must not be 
too large compared to unity. When amq ^ 1, the effects of the heavy quark doublers 
would need to be included in the chiral theory, and the analysis would become prohibitively 
complicated. A detailed discussion of this and other issues involved in incorporating heavy 



quarks into SXPT appears in Ref. 



19|. 



The calculations of interest here have been performed in continuum partially quenched 
chiral perturbation theory (PQXPT) by Becirevic, Prelovsek and Zupan 33| for Nsea de- 
generate sea quarks. In this paper we show how one can generalize the PQXPT formulae 
to the corresponding SXPT formulae, thereby avoiding the necessity of recomputing all the 
diagrams from scratch. 

Some results from the current work, as well as a brief discussion of how to generalize 
PQXPT to SXPT, appear in Ref. 1341. In addition, our results have already been used in 



chiral fits to lattice data in Refs. 



A related calculation for the B ^ D* and B ^ D 



semileptonic form factors has been presented by Laiho and Van de Water 28 1. 

The outline of this paper is as follows: We first include a brief description of heavy-light 
SXPT in Sec. [Tll In Sec. IIIIl we discuss the procedure for generalizing PQXPT to SXPT, 
using the heavy-light form factors as examples, although the procedure can be used for 



many other quantities in SXPT. Using this procedure and starting from Ref. [33|], we write 
down, in Sec. llVt the one-loop SXPT results for the semileptonic form factors. The partially 
quenched staggered case with non-degenerate sea quarks, as well as its continuum limit, is 
presented in Sec. IIV A[ In that section, we also discuss a method for treating — in a way 
that appears to be practical for fitting lattice data — some spurious singularities which 
arise in the calculations. Section IIVBI considers fuU-QCD special cases of the results from 



^ Taste violations with improved staggered fermions go like Q;|.a^. See Fig. 1 in Ref. 34] for a test of this 
relation. 



4 



Sec. IIV Aj while Sec. IIV CI discusses the analytic contributions to the form factors at this 
order. In Sec. |V] we add in the effects of a finite spatial lattice volume. Sec. |VT] presents 
our conclusions. We include three appendices: Appendix |X] gives expressions for the SXPT 
propagators and vertices, as well as the corresponding continuum versions. Appendix [B] 
lists the integrals used in the form factor calculations; while Appendix collects necessary 
wavefunction renormalization factors that were calculated in Refs. jisl . 3]. 



II. HEAVY-LIGHT STAGGERED CHIRAL PERTURBATION THEORY 

References 
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18j show how to incorporate heavy-light mesons into continuum 
XPT; the extension to SXPT appears in Ref. (lol. Here we review the key features needed 
for our calculations. 

The heavy-light vector (B^^) and pseudoscalar (Bo) mesons are combined in the field 

^« = ^ b^B;, + *755a] , (1) 

which destroys a meson. Here v is the meson velocity, and a is the "flavor-taste" index of 
the light quark in the meson. For n flavors of light quarks, a can take on 4n values. Later, 
we will write a as separate flavor (x) and taste (a) indices, a (x, a), and ultimately drop 
the taste index, since the quantities we calculate will have trivial dependence on the light 
quark taste. The conjugate field Ha creates mesons: 

Ha = loHho = [l,Bl: + t^.Bl] . (2) 

As mentioned in the introduction, we use B to denote generic heavy-light mesons to empha- 
size that we are working to leading order in l/mq. 

Under SU (2) heavy-quark spin symmetry, the heavy-light field transforms as 

H SH , 

H HS^ , (3) 
with S G SU{2), while under the SU^An)^ x SU{An)R chiral symmetry, 

H HU^ , 

H ^ VH , (4) 



with U G SU{An) defined below. We keep the flavor and taste indices implicit here. 

The light mesons are combined in a Hermitian field For n staggered flavors, $ is a 

4n X 4n matrix given by: 



u 


7r+ 


K+ 




D 




K- 




S 



\ 



(5) 



We show the n = 3 portion of $ explicitly. Each entry in Eq. ([5]) is a 4 x 4 matrix, written 
in terms of the 16 Hermitian taste generators T= as, for example, U = U~,T^. The 

component fields of the flavor-neutral elements (f/s, 1^=, . . . ) are real; the other (flavor- 
charged) flelds (vrj, K^, . . . ) are complex. The T= are 



(6) 



with C,fi the taste matrices corresponding to the Dirac gamma matrices, and = I the 4x4 
identity matrix. We deflne ^^,5 = ^^^5, and ^^.^ = (1/2) [^^, ^1,]. 
The mass matrix is the 4n x 4n matrix 

/ 



M 



mj 
mJ 
mJ 



V 



(7) 



where the portion shown is again for the n = 3 case. 

From $ one constructs the unitary chiral fleld S = exp [«$//], with / the tree-level pion 
decay constant. In our normalization, / ~ = 131 MeV. Terms involving the heavy-lights 
are conveniently written using use cr = = exp[i$/2/]. These flelds transform trivially 
under the SU{2) spin symmetry, while under SU{An)L x SU{An)fi we have 



(8) 
(9) 



with global transformations L G SU{An)L and R G SU{An)ji. The transformation U, deflned 
by Eq. ([9]), is is a function of $ and therefore of the coordinates. 
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It is convenient to define objects involving the a field that transform only with U and 
The two possibilities with a single derivative are 

= ^ [a^d^a + ad^a^] , (10) 
K = l Wd,a - ad,a^] . (11) 

transforms like a vector field under the SU{4n)L x SU{An)fi chiral symmetry and, when 
combined with the derivative, can form a covariant derivative acting on the heavy-light field 
or its conjugate: 

{D;H), = DfH, = d,Ha - i^fHb , (12) 

with implicit sums over repeated indices. The covariant derivatives and transform under 
the chiral symmetry as 

A^ ^ UA^U^ . (13) 

The combined symmetry group of the theory includes Euclidean rotations (or Lorentz 
symmetry), translations, heavy-quark spin, flavor-taste chiral symmetries, and the discrete 
symmetries C, P, and T. Many of these symmetries are violated by lattice artifacts and/or 
light quark masses. Violations to a given order are encoded as spurious in the Symanzik 
action. From these spurious, the heavy-light and light-light fields, derivatives, the heavy 
quark 4-velocity f^, and the light quark gamma matrix 7^, we can construct the chiral 
Lagrangian and relevant currents order by order. 

Reference 19| finds the lowest order heavy-chiral Lagrangian and left-handed current, 
as well as higher order corrections. We need primarily the lowest order results here. For 
convenience, we write the Lagrangian in Minkowski space, so that we can make contact with 
the continuum literature. 

We write the leading order (LO) chiral Lagrangian as 

^LO = -^pion + -^HL (14) 
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where £pion is the standard SXPT Lagrangian pj| for the hght-hght mesons, and £hl is the 
contribution of the heavy-hghts. We have^ 

£pion = ^IV(a^Ea^Et) + l/ifTr(>[S + A^St) 
o 4 

-'^-^{U, + D, + S, + ...f-a'V, (15) 
-V = CiTr(e5")Eef St)+C3i5^[Tr(e^Ee^E) + /..c.] 

+ C5V^ E T^e^^^S) Tr(e^Et) + cj- J2 Tr(ei?S) Trl^^^^Et) , (16) 

£hl = -I TiiHHvD) + Tt(HHY15K) ■ (17) 

Here Tr denotes a trace over flavor-taste indices and, where relevant, Dirac indices. The 
product HH is treated as a matrix in flavor-taste space: {HH)ab = HaHb. The covariant 
derivative D acts only on the field immediately preceding it. For convenience, we work with 
diagonal fields {U, D, . . . ) and leave the anomaly (wIq) term explicit in Eq. f|T5l) . We can 
take ml oo and go to the physical basis (7r°, f], . . .) at the end of the calculation js^]- 

To calculate semileptonic form factors, we need the chiral representative of the left-handed 
current which destroys a heavy-light meson of flavor-taste b. At LO this takes the form 

J^i = ^MY{l-j,)H)a^X^'K (18) 

where A*^''^ is a constant vector that fixes the flavor-taste: {\^''^)c = ^bc, and tr^, is a trace on 
Dirac indices only. 

The power counting is a little complicated in the heavy-light case, since many scales are 
available. Let rriq be a generic light quark mass, and let oc be the corresponding 
"pion" mass, with p its 4-momentum. Further, take k as the heavy-light meson's residual 
momentum. Then our power counting assumes /c^ ~ ~ ~ nig ~ a^, where appropriate 
powers of the chiral scale or Aqcd are implicit. The leading heavy-light chiral Lagrangian 



There is a missing minus sign in Eq. (35) of Ref. 19|. 
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Chl is 0{k), the leading light-light Lagrangian £pion is OijP', niq, a^), and the leading heavy- 
light current J^q is 0{1). Only these leading terms are relevant to the calculation of non- 
analytic "chiral logarithms" at first non-trivial order, which give 0{mg, a^) corrections to 
leading expressions for semileptonic form factors. 

In principle, finding the corresponding analytic corrections requires complete knowledge 
of the next-order terms in the Lagrangian and current. However, since the form factors 
depend only on the the valence and sea quark masses, a^, and the pion energy in the rest 
frame of the B (namely v-p), the form of these corrections is rather simple and is easily 
determined by the symmetries. The large number of chiral parameters that can appear in 
higher-order terms in the Lagrangian and the current collapse down into relatively few free 
parameters in the form factors. Unless one wants to write these free parameters in terms 
of the chiral parameters, complete knowledge of the higher-order terms in the Lagrangian 
and current is often unnecessary. However, one does need to know enough about the higher- 
order terms to check for the possibility of relations among the free parameters that multiply 
different quantities or that appear in different form factors. At the order we work here, there 
is one relation among the various parameters that determine the linear dependence of the 
two form factors on the valence masses. In order to be sure that this relation is valid, we 
need to know all terms at next order that can contribute such linear dependence. 

Fortunately, all such terms are known. For the light-light Lagrangian, Eq. f|T5l) . the 
relevant terms are the standard ~ m^) terms in the continuum jssj]. All terms of 

(9(mqa^,a^), which are special to SXPT, are also available 13]. For the heavy-light La- 

ler order than Eqs. ( ITTl) and (ITSl) 



19| does not attempt a complete 



grangian and current, Ref. 19|] lists all terms which are higl 
by a factor of nig (most important here) or a^. Reference 
catalog of the terms which are higher than Eqs. f|T71) and f|T8l) by one or two powers of k, 
i.e., having one or two derivative insertions. However, a sufficient number of representative 
terms of this type are listed to see that the corresponding free parameters in the form factors 
are all independent. We discuss the determination of the analytic terms further in Sec. llVCi 

III. GENERALIZING CONTINUUM PQXPT TO SXPT 

We wish to compute the decay —>■ Pr^y in SXPT, where x and y are (light) flavor 
labels. The taste of the light quarks in B, P and the current also needs to be specified. 
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We take the Pxy to be a "Goldstone pion" with taste ^5. Let the hght quark in the B 
have taste a {a = 1,...,4); in flavor-taste notation the hght quark has index a xa. 
The current, Eq. (IT5]) . has flavor-taste b ^ yf3. Despite the existence of taste violations 
at non-zero lattice spacing, the amplitude turns out to be proportional to (^5/2)^^, with a 
proportionality factor that is independent of the tastes a,f3. We will often keep this rather 
trivial taste-dependence implicit. 



In Ref. 331], Becirevic, et al. have calculated the form factors for B ^ tt and B —>■ K 
transitions in continuum PQXPT. They assume degenerate sea-quark masses, but leave Ns^a,, 
the number of sea quarks, arbitrary. As we explain below, the iVsea dependence is a marker 



for the underlying quark flow [29|] within the meson diagrams. Once we have separated the 
meson diagrams into their contributions from various the quark flow diagrams, we can easily 
generalize the continuum PQXPT results to the staggered case, without actually having to 
calculate any SXPT diagrams. To check our method, however, we have also computed many 
of the diagrams directly in SXPT; the results agree. 

The key feature that makes possible the generalization of continuum PQXPT results 
to SXPT results is the taste-invariance of the leading-order Lagrangian for the heavy-light 
mesons |19j]. This means that the continuum vertices and propagators involving heavy- 
light mesons are trivially generalized to the staggered case: flavor indices (which can take 
Nsea values if they describe sea quarks) simply become flavor-taste indices (taking 4A^sea 
sea-quark values). In one- loop diagrams, taste violations arise only from the light meson 
("pion") propagators. Propagators and vertices for the staggered and continuum cases are 
listed Appendix El 



Looking at the expressions in Appendix B of Ref. [33], we see that there are two types 
of terms that can contribute to each diagram for B^ — > Pxy'- a term proportional to A^^sea, 
and a term proportional to l/N^^.. This is the same behavior that appears, for example, in 



light-light [351] or heavy-light [361] PQXPT decay constants. 

The term which is proportional to Asea comes solely from connected quark-level diagrams, 
an example of which is shown in Fig. [1] (where (a) is the meson-level diagram and (b) is 
the quark-level diagram).^ The appearance of the quark loop accounts for the factor of 



By definition, the pion propagator in Fig. [DJa) is connected; the version with a disconnected propagator 
is shown in Fig. [D^a). 



10 



A^sea- In detail, using Eqs. (lAlSP and (lAlGp . the loop integrand is proportional to the 
connected contraction X]j{'^«i'^i*'}conn' ^^^^^ index j is repeated because the heavy- 
light propagator conserves flavor. Equation (lAlSp then implies that the sum over j produces 
a factor of A^^sea when the sea quarks are degenerate. In the non-degenerate case, there is no 
factor of A^sea but simply a sum over the sea-quark flavor of the virtual valence-sea pion. 

In the staggered case, the internal heavy propagators, Eqs. (lAip and (lA2p . as well as the 
vertices coupling heavy-light mesons to pions [e.g., Eq. flA3l) ). preserve both flavor and taste. 
Therefore Fig. [1] is now simply proportional to Y.b{^ab^ba'} ^^^^ = Y. j fsi^ f^^ j i' c^'} conn^ 
where we have replaced the flavor-taste indices (a, 6, ... ) with separate flavor . . . ) and 
taste (a, ) indices. From Eq. (lASp . the loop integrand is then proportional to 

y , '^''f'-' . , (19) 

where the 6aa' factor shows that, despite the existence of taste violations, the loop preserves 
the taste of the light quark in the heavy-light meson and is independent of that taste. 

The overall factor of the SXPT diagram must be such as to reproduce the continuum 
result in the a — > limit. Since pions come in 16 tastes, the sum over pion tastes S in 
Eq. fll9p must come with a factor of 1/16 compared to the continuum expression. To see 
this explicitly, note first that there are two factors of 1/2 relative to the continuum coming 
from the vertices (compare Eqs. flA3l) and flA16P ). due to the non-standard normalization of 
the taste generators, Eq. ([6]). An additional factor of 1/4 comes from the SXPT procedure 
for taking into account the fourth root of the staggered determinant: This is a diagram with 
a single sea quark loop. 

Finally, we need to consider how such a diagram depends on the tastes a and /? of the 
heavy- light meson and the current. Since the taste indices flow trivially through the heavy- 
light lines and vertices, and, as we have seen, through the loops, the taste dependence is 
simply {^^/2)ai3, where the .^5 comes from the outgoing light meson. The factor of 1/2 is due 
to the normalization of the taste generators. 



The net result is that terms with factors of A^^sea in the continuum calculation of Ref. 
are converted to SXPT by the rule: 



33| 



N...Hml) - ^ 1^ E ^H.^) (20) 



16 
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where the sum over / is over the sea quark flavors, z is the valence flavor flowing through 
the loop (either x or y), tum is the common mass of the Nsea mesons made up of a z valence 
quark and the degenerate sea quarks, and is some function of the pion masses. (For 
heavy-light quantities, JF is often also a function of the pion energy in the heavy-light rest 
frame.) The masses of pions of various tastes and flavors (ttij^, =) are given in Eq. (lA7p . 

The terms that are proportional to 1/Nsea are more subtle. They arise from diagrams with 
disconnected pion propagators. The simplest example is shown at the meson level in Fig. [2t^a) 
and at the quark level in Fig. [2](b). The continuum form of the disconnected propagator is 
given in Eq. ( ]A21I) . Using the continuum values S' = mg/S and m'^ ^ NscaTnl/3, we see that 
the disconnected propagator produces an overall factor of 1/A^sea as rriQ oo. Equation 
f lA2ip can then be written as a sum of residues times poles, where the residues can be rather 
complicated when the sea masses are non-degenerate (see Appendix [B]) . Thus, the flnal 
answer after integration amounts to something of the form 



^-J^R^Hm',), (21) 



iVsea 

3 



where JF is again a general function resulting from the loop integral, Rj is the residue of 
the pole at = m|, and j ranges over the flavor-neutral mesons involved: the sea mesons, 
710,7],..., and the "external" mesons in the disconnected propagator, called ii and i'i' in 
Eq. (]A2ip . When mj/j/ = rrin, there is a double pole, and Eq. (!2T|) should be replaced by 

]^E5^[«.-^K)] . (22) 

where the sum over j now does not include i'i'. 

When the sea quarks are degenerate, the residues simplify considerably. However, by 



comparing the general forms in Eqs. ( 12T|) and ( l22l) to the rather simple terms in Ref. [33|], it 
is easy to move backwards from the degenerate case and determine the form of the expressions 
for non-degenerate sea quarks. 

The flavor structure in the staggered case is identical to that in the continuum: Flavor 
remains a good quantum number, so meson propagators in both cases can only be discon- 
nected if they are flavor neutral. Because of taste violations, however, disconnected hairpin 
diagrams can contribute to mesons propagators with three different tastes (singlet, vector, 
and axial vector) at this order in SXPT. These three hairpin contributions are quite similar 
to each other, but there are a few important differences: 

12 



• The strength of the hairpin, (54, depends on the taste S — see Eq. (lA9p . 

• In the taste-singlet case, as in the continuum, the hairpin {4171^/3) comes from the 
anomaly and makes the flavor-singlet meson heavy. Decoupling the rij by taking 
ttiq ^ oo is therefore a good approximation, and we do it throughout, giving rise to 
an overall factor of 1/A^sea- But in the taste- vector and taste-axial- vector cases, the 
hairpins are not particularly large; indeed they are taste-violating effects that vanish 
like (up to logarithms) as a ^ 0. So we cannot decouple the corresponding mesons, 
r]y and r/^, in the taste- vector and taste-axial- vector channels. 

• The taste matrices associated with the vector and axial- vector mesons, and ^^5, 
anticommute with the ^5 coming from the outgoing Goldstone pion. Therefore the 
vector and axial hairpin contributions will have an opposite sign from the singlet (and 
continuum) contribution if the ^5 needs to be pushed past a or ^^5 to contract with 
the external pion state. 

Figure [3] shows the tree-level diagrams that contribute to the form factors, while Figs. H] 
and [5] show all the non-vanishing one-loop diagrams. As a first example of the treatment 
of diagrams with disconnected meson propagators, consider Fig. [5](b). It is not hard to see 
that this diagram has only a disconnected contribution, shown as a quark-fiow diagram in 
Fig. El A connected contribution would require the contraction of the external light quark 
fields X and y, which make up the outgoing pion. That is impossible since we have chosen 
X 7^ y.^ 



In our notation, the result from Ref. 
quenched case with A^sea degenerate sea quarks is: 

9l 1 



for this diagram in the continuum partially 



"''^ jr (my , VP) - "^h rimx, v-p) 



2 2 1 V" "-T ; ^ fj 2 2 

my — mj^ ttly — rrix 



(23) 



(47r/)2Ar„ 

where mu is the mass of any of the mesons made up of a sea quark and a sea anti-quark, 
mx and my are the masses of the fiavor-neutral mesons made up of xx and yy quarks, 
respectively, and the function J™*^, defined below in Eq. ( l5Ti) . is the result of the momentum 
integral. 



A similar argument will be given in more detail below when discussing Figs. [7] and [51 
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The ratios of mass differences in Eq. (123|) can be recognized as the residue functions (see 
Appendix [B|) for the various poles. For example, (my — m^)/(my — m'j^) is the residue for 
the pole at = my. These residues are rather simple in this case because of the degeneracy 
of the sea quarks. To generalize Eq. fl2^ to the completely non-degenerate case, we simply 
need to replace the residues by their general expressions. For Ngea. non-degenerate sea quarks, 
Eq. (123|) is replaced by 

- j^it ^ [«f ■'r''(",. -p)] . (24) 

where the Minkowski residues ^^"'^^ are defined in Eq. ( ]B2[) . and the sum over j is over the 
Nsea, + 1 mesons that make up the denominator masses in the disconnected propagator after 
ml — > oo. (See Eq. (1A21I) and the discussion following it.) We leave implicit, for now, the 
arguments to the residues in Eq. (!24|) : we will be more explicit in the final results below. In 
addition, we will ultimately express everything in terms of Euclidean-space residues 
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19|. 



Eq. (]B4p . simply because those are what have been defined and used previously 

Cases with double poles present no additional problems, since Ref. [331] shows these 

explicitly as derivatives with respect to squared masses of the results of single-pole integrals. 

We will therefore simply get derivatives of the usual residues, as in Eq. ( ]B3[) . 

As discussed above, we will need the expression before the itlq oo limit is taken in order 

to generalize the result to the disconnected taste-vector and axial-vector cases. Equation 

(IA21|) and the fact that m'i, ^ NscaiTil/S for large mo allow us to rewrite Eq. flM|) as 



+(if|pf 5:[^r"''"-''-^r''(™„-p)] . (25) 

The sum over j now includes the rj'. The sign difference between Eqs. f l24p and (123]) comes 
from the sign of the mass term in the Minkowski-space r]' propagator. 

Generalizing Eq. to the staggered case is then straightforward. For the taste-singlet 
hairpin contributions, we simply replace each continuum pion mass by the mass of the 
corresponding taste-singlet pion. In other words, we just let rrij —>■ rrijj in Eq. fl25|) . Note 
that, after the staggered fourth root is properly taken into account, the taste-singlet rj' mass 
goes like for large mo, as it does in the continuum, so one could reverse the process 

that led to Eq. (123]) and use instead Eq. or even Eq. (12^ (for degenerate sea-quarks), 
with rrij — > rrijj in both cases. Just as for diagrams with connected pion propagators (see 
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Eq. (120|) ). there is also a trivial overall factor of {^5)013/2, where a and /3 are the tastes of 
the heavy-light meson and the current, respectively, and the ^5 is due to the pseudoscalar 
(Goldstone) taste of the outgoing pion. 

For the taste-vector and axial-vector disconnected contributions, a little more work is 
required. We first note that the factor of 171^/3 in Eq. fl2^ is simply S'^/A with 3 = 1, 
the strength of the taste-singlet hairpin, Eq. (lA9p .^ For the other tastes we then replace 
54; by the appropriate hairpin strength from Eq. flA9|) and also replace the pion masses: 
rrij ^j,-.- In addition, there is an overall sign change for this diagram in going from the 
singlet to the vector or axial-vector tastes. This comes from the fact that the outgoing pion 
line in Fig. El^b) lies between the two ends of the disconnected propagator. Using Eq. (1A13P 
and the Feynman rules for the heavy-light propagators and vertices in Appendix |Xl one sees 
that the diagram with a taste-S disconnected propagator goes like ( Th ^5 Th )^^. This leads 
to a positive sign for S = J but a negative sign for tastes that anticommute with ^5. Finally, 
the fact that there are four degenerate taste- vector (or axial-vector) pions at this order leads 
to an additional overall factor of four. 

When we attempt to apply the same procedure to the other diagrams in Figs. H] and [5l 
we find a further complication in diagrams Fig. Hl^a), Fig. IH^b), and Fig. [5]^c), where the 
external pion and one or more internal pions emerge from the same vertex. The problem is 
that the ordering of the taste matrices at the vertex is not determined by the meson-level 
diagram {i.e., each diagram can correspond to several orderings), so we do not immediately 
know the relative sign of taste-vector and axial-vector contributions relative to the singlet 
contribution. Nevertheless, a quark-fiow analysis allows us to identify appropriate "fiags" 



that signal which terms in Ref. j33j come from which orderings at the vertex. 

As an example of the procedure in this case, consider Fig. [5](c). The corresponding 
quark fiow diagrams with disconnected pion propagators are shown Fig. [71 In Fig. [7](a), 
the outgoing pion lies between the two ends of the disconnected propagator. This produces 
a change in sign of the taste-vector and axial-vector hairpin contributions relative to the 
taste-singlet one, just as for Fig. [5](b). In Fig. [7](b), on the other hand, the outgoing pion 
is emitted outside the disconnected propagator, and all the hairpin contributions have the 



^ The factor of 1/4 just comes from the different conventional normahzation of the generators in the con- 
tinuum and staggered cases; see Appendix VK\ for further discussion of normahzation issues. 
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same sign. The same is true of the reflected version of Fig. [7](b), which has the outgoing 
pion emerging from the other side of the vertex. 

Fortunately, Figs. [7](a) and (b) are distinguished by their flavor structure, even in the 
continuum. In Fig. Wi^), the two "external" mesons in the disconnected propagator have 
different flavors: The one on the left is an X meson (an xx bound state); while the one on the 
right is a y meson (a yy bound state). In Fig. Wljo), both external mesons in the disconnected 
propagator are Y mesons. Similarly, the reflected version of Fig. [7|(b) has two X mesons in 
the disconnected propagator. This flavor structure is immediately apparent in the results of 
Ref. [sS]. The parts of Fig. [5](c) that come from the quark flow of Fig. [Tl^a) are proportional 
to the function called Hi, which depends on the masses mx and my (in our notation), 
as well as the sea-meson mass. The parts of Fig. O^c) that come from the quark flow of 
Fig. Et^b) (or its reflected version) are proportional to the function called Gi, which depends 



only on the mass my (or mx) and the sea- meson mass. To generalize the results of Ref. [33 1 
to the staggered case, we thus can use the method outlined above, and simply include an 
extra minus sign for those taste-vector and axial-vector hairpin contributions proportional 
to Hi (relative to the taste-singlet contributions), but not for those proportional to Gi. This 
approach also works for the other problematic diagrams. Figs. 111(a) and (b). 

The reader may wonder why the complication associated with ordering the taste matrices 
at the vertices does not occur when the internal pion propagator is connected, but only in 
the disconnected case. Figure M shows possible quark-flow diagrams for Fig. M^c) with a 
connected pion propagator. Figure (Hl^a) cannot occur in our case because we have assumed 
that X, the light flavor of the heavy-light meson, is different from y, the light flavor of the 
weak current.''' The same reasoning is what allows us to rule out any connected contributions 
to Fig. [5t^b), as mentioned above. Thus all contributions with connected propagators are of 
the type shown in Fig. [8](b), or its reflected version, and these never have a sign difference 
between terms with different internal pion tastes. 

We note that one can reproduce the SXPT results for light-light llSj] and heavy-light 
jigj l mesons by starting from the continuum PQXPT in Refs. 35 1 or j36[, respectively, and 

following the procedure described above. The computations are in fact slightly more difficult 

ii 

in those cases than in the one at hand, because Refs. [35] and [36] do not explicitly separate 



^ Equivalently, we have assumed that the outgoing pion is flavor charged. 
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double-pole from single-pole contributions. It is therefore takes a little work to express the 
answers from those references in the form of our residue functions, which is the necessary 
first step before generalizing to the staggered case. 



IV. FORM FACTORS FOR B ^ P DECAY 



The standard form factor decomposition for the matrix element between a meson and 
a Pxy meson is 



iPB+p)t, - g^- 



rrin — m 

Ox 



2 



g2 



F^{q') + ^^-^q,F,{q') , (26) 
Q 



where q = ps — p is the momentum transfer. We are suppressing taste indices everywhere, 
but emphasize that the light pseudoscalar P^y is assumed to be the Goldstone meson (taste 
^5). In the heavy quark limit, it is more convenient to write this in terms of form factors 
which are independent of the heavy meson mass 

{Pxy{p)\yj,,b\Bx{v))j^Qj^j, = [Pm - {v-p)v^,] fp{v-p) + vjy{v-p) , (27) 

where v is the four- velocity of the heavy quark, and v-p is the energy of the pion in the 
heavy meson rest frame. Recall that the QCD heavy meson state and the HQET heavy 
meson state are related by 

\B{pb))qcd = \/m^\B{v))HQET ■ (28) 

The form factors fp and and are often called f± and /y, respectively. As discussed in 
Sec. the taste indices are left implicit in Eqs. f l26l) and ( 127|) . as are the trivial overall 
factors of (^5/2) in the matrix elements. 

The tree-level diagrams for B^ Pxy are shown in Fig. [31 Figure [3](a) is the tree-level 
"point" contribution to f^, while Fig. [3](b) is the tree-level "pole" contribution to fp. We 
have 

/rM = y, /r(M = 7^;:^' (29) 

where A* = itlb* — ms is the mass difference of the vector and pseudoscalar heavy-light 
meson masses at leading order in the chiral expansion, i.e., neglecting all effects of light- 
quark masses. As in Refs. 0, we drop this splitting inside loops, but keep it in the 
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internal B* line in the tree-level diagram Fig. [3](b). This forces the tree-level pole in fp to be 
at rriB*, the physical point. Dropping A* inside loops is consistent at leading order in HQET, 
which is the order to which we are working. It would also be consistent, parametrically, to 
drop A* everywhere. But this would not be convenient, since the m^. pole is physically 
important for fp. 

The non-zero diagrams that correct the form factors to one loop are shown in Fig. H] for 
fy and Fig. [5] for fp. Table [T] lists the correspondences between these diagrams and those 



of Ref. 



33| . A number of other diagrams, which can arise in principle, vanish identically 



due to the transverse nature of the B* propagator, Eq. (jA2]); these additional diagrams can 



be found in Ref. 



331] . We do not indicate hairpin vertices explicitly in Figs. H] and the 



internal pion propagators in these diagrams may be either connected or disconnected. 
Before generalizing the results in Ref. [33| to SXPT, we discuss a subtle issue that affects 



Fig. [5]^a). If the splitting A* is dropped on internal B* lines in loop diagrams, as is done 
in Ref. t33|, this diagram has a spurious singularity (a double pole) at v ■ p = 0, the edge 
of the physical region. The singularity arises from the presence of the two B* lines that 
are not inside the loop integral and therefore can be on mass shell in the absence of B*-B 
splitting. Including A* on all such internal "on-shell" B* lines {i.e., lines not inside the loops 



themselves), as is done in Ref. 



40|, at least pushes the unnatural double pole out of the 



physical region. We will follow this prescription for including the splitting, but take it one 
step further. The loop in Fig. [5]^a) is a self-energy correction on the internal B* line. The 
double pole results from not iterating the self-energy and summing the geometric series. 
We will follow the more natural course and sum the series; doing so restores a standard 
single-pole singularity. 

There is a further one- loop contribution that can naturally be included in Fig. [5]^a). The 
corresponding tree- level graph. Fig. [3]^b), gets two kinds of corrections that are not shown in 
Fig. [5l One comes simply from the wavefunction renormalizations on the external pion and 
B lines; we include those terms explicitly below. The second contribution arises from the 
one-loop shift in the external meson mass. Since this mass shift depends on the flavor of the 
light quark in the B^, namely x, we call it 5Mj. = H^iv ■ k = 0) , with ^^^.(f ■ k) the self-energy 
for B^ or B*. Note that Sj, is the same for both the B^. and the B*, since the splitting 
A* is dropped inside loops. When the external B^ line in Fig. [3](b) is put on mass-shell 
at one loop, the denominator of the internal B* propagator changes from —2{v-p + A*) to 
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-2{v-p + A* — SMx). It is convenient to define A*^ as the full splitting between a B* and a 



il, = Mb^ - Mb^ = a* + 5My - 



(30) 



The internal B* propagator now becomes —2{v-p + A*^ — 6My). The contribution from the 
mass shift may then be combined with the tree-level and Fig. El^a) contributions to give: 



K 



9-K 



D{vp) = S,(t;-p)-S,(0) , 



(31) 
(32) 



where the subtraction in D comes from the effect of putting the B^ on mass shell, via 
Eq. (EOD. 

The main difference between the approach taken to the spurious singularity of Fig. [5t^a) 



and that of Becirevic et al. [33|] is that they work to first order in the self-energy in the 
corresponding diagram (their diagra m ( 7)). Expanding Eq. ( !3TI) . we find that D is related 



in the continuum limit to what Ref. [33[ calls 5fp by 



D{vp) = -VP Sff^ . 



(33) 



(7) 



Thus we can find the staggered D{v-p) simply by applying the methods of Sec. Illll to 6f] 

We can now write down the expressions for the form factors for B^ — » P^y decay. For the 
point form factor, we have 



/^^ [1 + + c>,. + c>, + cL(m„ + m, + m,) 

+ divp) + clivp)' + cy] , 



(34) 



where fl^^'^'^ is given by Eq. ( !29l) . and the analytic coefficients from next-to- 

leading order (NLO) terms in the heavy-light chiral Lagrangian (see Sec. IIVCI) . The non- 
analytic pieces, which come from the diagrams shown in Fig. H] as well as the wavefunction 



renormalizations, are included in 5fv 



SfB.^p.. = + + 16Zb. + \5Zb., 



(35) 



The wavefunction renormalization terms, 5Zb^ and SZp^ , have been calculated previously 



13|, 



19| in SXPT and are listed in Appendix O 
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For the fp form factor, we write 

+ c?(t;-p) + c|(t;-p)2 + cy] . (36) 
where fp^^^ is defined in Eq. f lHTi) . and 

^"=<"-' - 7 ^ ■ 

Non-analytic contributions are summarized in the function D{v-p) in fp'^^^, Eq. (152]) . and 
5fp , which comes from Figs. [5]^b)-(d) and wavefunction renormahzations. Exphcitly, 

For simphcity, we do not include the superscript P^y on the individual diagrams in 

Eqs. ([35]) and dag). 

Using ff^", which includes the full splitting A*^, rather than Eq. (|29l), 

changes Eq. ( !36l) only by higher-order terms. However, it is convenient to keep the same 
splitting in both fp^^^ and the other terms in Eq. (|36ll . Note that it is also consistent at this 
order to use the alternative form 

fBx^Pxy ^ fseM ^ 5fB.^P^y + ^^^^ + + c^,,(m„ + + m.) 

+ c?M + c|(t;-pf + c^a^] , (39) 

The analytic terms in and fp are not all independent. As mentioned in Sec. [Tll there 
is one relation among the terms that control the valence mass dependence: 

c^ + c^ = c^ + c^ (40) 

We show that this relation follows from the higher order terms in the Lagrangian and current 
in Sec. IIVCI All other NLO parameters in Eqs. (134|) and (l36l) are independent. 



A. Form factors for 3-flavor partially quenched SXPT 

First we display the results for the individual diagrams shown in Figs. H] and [5] for the 
fully non-degenerate case with three dynamical flavors (the "1 + 1 + 1" case). This means 
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that we have aheady taken into account the transition from 4 to 1 tastes per flavor. Indeed, 
our method of generahzing the partially quenched continuum expressions to the staggered 
case automatically includes this adjustment. We detail below the minor changes needed to 
obtain 2 + 1 results from those in the 1 + 1 + 1 case. 

We first define sets of masses which appear in the numerators and denominators of the 
disconnected propagators with taste labels implicit (see Appendix [B|) : 



_y\/J (5,a;j/) 



{ml, ml, ml} 



r 2 2 2t 

{mx,m^o,m^} , 



2 _2 _2i 
X 

r 2 2 2 2 1 

{mx,m^o,m^,m^,\ , 

r 2 2 2 21 

{mx,mY,m^o,m^} , 

r 2 2 2 2 2 1 

{mx,mY,m^o,m^,m^,} 



(41) 
(42) 
(43) 
(44) 
(45) 



For the mass sets (H2l) and ( H3l) . there are also corresponding sets with x ^ y and X ^ Y. 
When we show explicit taste subscripts such as / or on the mass sets fi or Ai, it means 
that all the masses in the set have that taste. 

The functions that appear in the form factors are^ 

,2- 



Ji (m) = In ( — — 



m' 



Urn, A) = -2A^ In - 4A^F (^) + 2A^ , 

Ji(m,A) = (-m2 + ^A2)ln(^)+^(A2-m2)F 

Vl -a;2 tanh"^ (^Vl - x"^^ , < a; < 1 
— Va;^ — 1 tan~^ x^ — 1 j , x > 1 . 



10 A 2 4 2 
— A^ + -m^ 

9 3 



F{x) = 



(46) 
(47) 

(48) 

(49) 



The main difference in these formulae with those of Becirevic et a l. |33l| is that they keep 



the divergence pieces, while we have renormalized as in Refs. [13|, Il9|. To convert to our 
form, replace the MS scale in Ref. 33|] with the chiral scale A and set their quantity A to 
zero, where 



A 



4-d 



7 + ln(47r) + 1 



(50) 



For ease of comparison to Ref. 33|, we use /i(m) instead of £(rn^) (as in Refs. [3, [3]) for the chiral 
logarithm. 
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with d the number of dimensions. F{x) is only needed for positive x; so we use the simpler 



40|, rather than the more general version worked out in Ref. 4l|] and 



form given in Ref. 

quoted in Ref. We do not list the function J2, which appears in the integral J^" of 

Eq. fIBSp but does not enter the final answers. 
We also define a "subtracted" Ji function by 

27rm^ 



jr'(m,A) = Ji(m,A) 



3A 



(51) 



The subtraction term cancels the singularity when A ^ 0. The function Jf^^ enters naturally 
in the expression for the self energy correction D{v-p) because of the the subtraction in 
Eq. fl32l) . It also turns out to arise from the integral in Fig. [5](b) — see Eq. (26) in Ref. 40|. 
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For the point corrections in the 1 + 1 + 1 case, we have 



1+1+1 



1+1+1 



+ 



J2 Rf'^ (-M^^^^; /if ) [him,,j) + 2/2(m,,„ vp)] 



d 



dm 



YJ 



+ a^5 



2x1 
V 



d ( 



dmyy 



E i?f^(-M?'^^/if)[/iKy) + 2/: 



2[mjy, 



- 5^ ^1 (A^f /if) [/i(m,y) + 2h{m,y, vp)] 

j£j^{5,xy) 

+ [V-.A] 

' ( E i?f^^(A^f^^/if)/iK-.)) 



+ 



dmyj 



• j^_\4(3,y) 



( Y. fl;"'(A^^^''';/.S?')/■("J 

\ -.-^ A 7/1 



jg_A4(4,a:y) 



+ a'6 



2x' 



dmyy 



+ 



9 / 



( X: flr'(Air';M?)A('%v)) 

+ |V A| 



23 



Those that correct the pole form factors are 

3glv-p J 1 



l-^ J 1+1+1 



1+1+1 



^ / 1+1+1 



1+1+1 



+ 



(47r/)2 ] 16 
1 



d 



dm 



YJ 



R 



[3,3] 



+ a% Yl 



d 



dm 



R 



4,3] 



9l 



6(47r/)2 I 16 

d 



+ 



jgX(3.!/) 



9 



(3,y). ,,(3) 



R 



[3,3] 



+ 2 [R^;^'^{Mr^\,r)l,im,,) 

9 



I 2 r/ 

+ a 



xy 



(4,x). „(3) 



2 5: [flfl(Air»';f,i?')/,(m,v.) 



+ [V^A] 



9 



3 ^ dmt T 



+aH'y J2 



d 



dmyy 



(54) 



+a25t. 5^ ^] [Mf'''^■,^^f) Jr\m,y,vp) + [V A] 1,(55) 



(56) 



[V-A] .(57) 
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In Eqs. (152|) through (157|) . the exphcit factors of 1/3 in front of terms involving the taste 



singlet (/) mesons come from the factors of 1/Nsea in Ref. 33 1. 

To get the full corrections for both and /p, we need to add in the wavefunction renor- 
malizations, given in Appendix O in Eqs. fl5^ and fl55]) . Putting these together with the 
analytic terms and (for fp) the D term, Eqs. (1341) and (I36p give the complete NLO expressions 
for the form factors in SXPT. 

The above 1+1+1 results are expressed in terms of the Euclidean residue functions R^p''\ 
Eq. (1B4I) . In the 2 + 1 case, there is a cancellation in the residues between the contribution 
of the U OT D in the numerator and that of the ttq in the denominator. Thus, to obtain the 
2 + 1 from the 1 + 1 + 1 case, one must simply reduce by one all superscripts on the residues, 
i.e., /^[""i./s-i]^ and remove and (say) uid from the mass sets: 

{m^,m|} , (58) 

A<(3,x) ^ {m^,mj}, (59) 

^^(4,.) ^ {ml,ml,m'^,} , (60) 

^ {m^,m^,mj} , (61) 

^^(^'"2^) ^ {m^,m^,mj,mj,} . (62) 

We also write here the expressions for three non-degenerate dynamical flavors in contin- 
uum PQXPT, which to our knowledge do not appear in the literature. These expressions 
can be obtained either by returning to Ref. 3^ and using the residue functions to generalize 
to the non-degenerate case, or simply by taking the continuum limit of the above equations. 
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Either way, the results for are 



cont 



2{Anf) 
1 



1 

+ 3 



9 



dmy 
d ( 



jgX(3,!/) 
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while those for fp are 

^^cont^ Bx^F 



cont 



cont 



cont 



f 



{Anf y ] 3 
1 



dmy 



6(47r/)2 |X^ [/i(m,.^) + /i(m,;)] 



1 

+ 3 



F — 

9 



Y 



[fiM(A1<''''';/'l)^i(m, 



— dm'ir 



+ 2 Y 

1 

^2(477)^ 



rrin 



-5 E a^K'"(^'''"'^^''')^.' 



d 



rrii 



(64) 



jgX(3,H) 

Corresponding continuum-hmit results for the wave-function renormalizations are given in 
Appendix O 



B. Full QCD Results 

Adding together the complete results for the "full QCD" case is straightforward. For 
simplicity, we specialize to case m„ = {i.e., 2 + 1). For B ^ it, the complete corrections 
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(including wave-function renormalization) are: 

je{n,'n,ri'} 

+ [V^A]\, (65) 



(47r/)2 I 16 4^ 



1 + 3^^ 



[2/i(m^,s)+7i(m;f,s)] 

1 1 1-1- Sn^ 

--^^Jr'^(m.,„^-p) + -^^Jr^(m,,„^;.p) + 

je{7r,»7,r)'} 



(66) 



(47r/)2 1 16 5 



[2/i(m^,s) + /i(mi^,s) 



+ 2/2(m^,s, ^;-p) + l2{mK,E, v-p) 



+ 



1 + 3^^ 



+ E 

jG{7r,r/,T;'} 

3(^^ - 1) 



a^6yRf'^^ {{m^y, m^,v, m^^y}; {msy}) 



X 



h{mjy) - 2l2{mjy,v-p) 



+ [V^A] 



(67) 
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For B K,^ we have 



.2 N d[3,1] 



+ ^ {-a'^Sy)Rf^^ {{msy, m^y, m^/y }; {m^y}) Jl''^{mjy, vp) 



j&{s,v,v'} 



+ [V-^A]\, 



(68) 



5/; 



B^K 



-y 

(47r/)2 I 16 4^ 



je{n,ri,r]'} 
je{5,7y,r)'} 



9l 



2 2 



(69) 



^ The transition B ^ K occurs through penguin diagrams; D K is a standard semileptonic decay due 
to the current in Eq. . We keep the notation B ^ K however to stress that we are working to lowest 
order in the heavy quark mass. 
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B^K 



(47r/)2 I 16 4^ 



+ 2l2{mK,E, v-p) + h{ms,r., v-p) 



12 



h{mr,j) + l2{mr,j,v-p) - l2{msj,v-p) 



V 



h{mrj>y) - h{m^,v) + h{mr,'y, v-p) - hirriny, v-p) 



" ^ Rf'^H{'i^sy,mr^y,mr,'y}]{m^y}) (^himjy) + l2{mjy,v-p)j 

+ 5Z ^5^'^' ({"^TT.V, "^ryy, "^r,',y}; {"^5,^}) hinijy) 



+ [V^A] 



(70) 



C. Analytic terms 

From the power counting discussed in Sec. [Tll as well as interchange symmetry among the 
sea quark masses, the form factors at the order we are working can only depend only on the 
valence quark masses rrix and rriy, the sum of the sea quark masses rriu + Trid + ms, the pion 
momentum (through v-p), and the lattice spacing, a. The last must appear quadratically, 
since the errors of the staggered action are 0{a'^). Recall that we do not include any 
discretization errors coming from the heavy quark in our effective theory. 

Thus we expect to have the analytic terms shown in Eqs. flM|) and fl36l) with coefficients 
cf and c^. (Here i = {x, y, sea, 1, 2, a}.) We then can examine, one by one, the known NLO 
terms in the Lagrangian and current to check for the existence of relations among the 
and/or c^. As soon as a sufficient number of terms are checked to ensure that the parameters 
are independent, we are done. It is therefore not necessary in all cases to have a complete 



catalog of NLO 
come from Ref. 



terms. Unless otherwise indicated, all NLO terms discussed in this section 



19|. 



Note first of all that we do not need to include explicitly the effects of mass- 
renormalization terms in the NLO heavy-light Lagrangian, such as 



2Ai Tr {HHM+) + 2\[ Tr (HH) Tr (M^) 



(71) 
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where we define 



A1± = 1 [aMa ± a^Ma^) . (72) 



The effect of the terms in Eq. ( ITTj) is absorbed into the B*-Bx mass difference A*^, Eq. ( !30l) . 
just hke the one-loop contribution to the mass. Corresponding 0{a?) term in the Lagrangian, 
which can be obtained by replacing M.^ above by various taste- violating operators, can 
likewise be ignored here. 

We now consider the discretization corrections parametrized by and c^. There are a 
large number of 0{a^) corrections to the Lagrangian and the current that can contribute to 
these coefficients, so it is not surprising that they are independent. For example, consider 
the following terms in the NLO heavy-light Lagrangian 



a" 



5^c^,Tr(M7,75{A^O,^•+}) , (73) 
fc=i 

where the O^'^ are various taste-violating operators, similar to those in Eq. (fT6l) above. 
These terms do not contribute to c^, but only to c^, though corrections to the B-B*-7r 
vertex in Fig. [3]^b). On the other hand, there are many terms that contribute both to 
and to c^. An example is the following correction to the current 

a' E tr, (7^(1 - l,)H) Ot'^-cr^X^'^ , (74) 

k=l 

which contributes equally to and c^. Additional examples are provided by those terms 



with two derivatives in the 0{mqa?) pion Lagrangian 39], which correct both coefficients 
though their effect on the pion wave-function renormalization. 

We consider the v-p and {v-pY terms next, namely c^, Cg, c^, and c^. This is a case 
where a complete catalog of Lagrangian and current corrections does not exist. However, it 
is easy to find corrections that contribute only to or only to fp. As in the previous case, 
corrections to the B-B*-^ vertex in Fig. [3](b) only affect fp at the order we are working. 
Thus, 

^ Tr ({v ■ DHH - HHv ■ D) 7^754^] (75) 



contributes to only; while 



Tr(m7^75(^7-:D)=^A^) (76) 



A2 

X 
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contributes to only. Similarly, only is affected, though Fig. [3](a), by any correction to 
the current whose expansion in terms of pion fields starts at linear order {i.e., corrections of 
schematic form + ■■■), with • • • denoting higher order terms in $). Thus, 

^tr,(7'^(l-75)i^)t;-AaU(^) (77) 

contributes to dl only; while 

^ tr,(7^ il-^,)H)v ■ Dv ■ AatA^") (78) 

X 

contributes to C2 only. Since there is at least one Lagrangian or current term that contributes 
to each of c^, Cg, c^, and exclusively, these coefficients are independent. 

The argument for the independence of the sea-quark mass terms, i.e., the coefficients c^^^ 
and cfg^, is similar. The Lagrangian correction 

h Tr (HH^^^.A^) Tt{M+) (79) 

contributes to cf^^ only; while the current correction 

P2 tio - 75)^) cr^X'^'^ Tr(A^+) (80) 

contributes equally to both cfg^ and Cgg^- These two observations are enough to guarantee 
that Cggg^ and c^^^ are independent. 

We now turn to the coefficients that control the valence quark mass dependence of the 
form factors: c^, Cy, c^, and c^. At first glance, it would seem unlikely that there could 
be any constraint among these parameters, since there are seven terms in the Lagrangian 
and current in Ref. |[19j] that could generate valence mass dependence. However, three of 
these terms are immediately eliminated, either because they could only contribute to flavor- 
neutral pions (with x = y), or because they produce no fewer than two pions. There are 
then two remaining corrections to the heavy-light Lagrangian, 

iki Tr (hHvDM^ - vDHHM^^ + h Tr (HH^^^r^{M, M^}) (81) 

and two corrections to the current, 

pi tr, (7^(1 - ^,)H) A^+ortA^") + p3 tr, (7^! - 1^)H) M-a^X^'^ (82) 

-"^^ There are additional terms involving Tr(A^ + ), as in Eqs. (I79p and ((80|) . that only give sea quark mass 
dependence at this order. 
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The fcs term in Eq. (IHT]) contributes only to fp, through the B-B*-tt vertex. However, 
because of the anticommutator, its contribution is proportional to nix + my, so it gives equal 
contributions to and c^. Similarly, because the one-pion term in is proportional to 
^M. + A^$, the ps term contributes equally to c!^ and (but not at all to and cp. 
Further, since M.^ creates only even number of pions, we can replace it by M. in Eq. (1821) . 
The pi term can then easily be seen to contribute equally to and c^, since the current 
needs to annihilate a B* in the fp case and a B^ in the case. 

The contributions of the ki term in Eq. (IHTj) are the most non-trivial. It contributes to 
both and through wave function renormalization on the external B^ line, but it also 
contributes to through an insertion on the internal B* line in Fig. [3](b). However, since 
wave-function renormalization effects on external lines go like a/Z, the contributions of this 
term to both and eg are exactly half of its contribution to c^. Thus, all four terms in 
Eqs. ( IHTj) and (!82|) are consistent with the relation given in Eq. (140|) . 

We still need to worry about valence mass dependence generated by the standard 0{j)^) 
pion Lagrangian 38|] through wave function renormalization of the external pion. Such 
contributions do exist (from L5), but the x ^ y symmetry of the pion guarantees they are 
proportional to m^. + ruy in both fp and /„, and hence do not violate Eq. fj40|) . 

A consistency check of the relation, Eq. (HOj) . as well as of the claimed independence of the 
other analytic terms, can be performed by considering the change in the chiral logarithms in 
Eqs. (|52ll through (!57ll and Eqs. ( IClll and ( 1C2I) under a change in chiral scale. To simplify 
the calculation, it is very convenient to use the conditions obeyed by sums of residues, which 
are given in Eq. (38) of the second paper in Ref. 13[]- We find that such a scale change can 
be absorbed by parameters that obey Eq. fHOj) but are otherwise independent. 

In the continuum limit, c^^^ and c^^^ remain independent, as do c^, c^, c\, and Cg. We 
disagree on these points with Ref. 3^, which found c^^^ = c^^^, and did not consider analytic 
terms giving v-p dependence. The difference can be traced to the inclusion here of the effects 
of the complete set of NLO mass-dependent terms, as well as a sufficient number of higher 
derivative terms (Eqs. ( 175|) through (1781)). In particular, the independence of c^^^ and c^^^ can 
be traced to the existence of the Lagrangian correction, Eq. f l79|) . which was not considered 
in Ref. 33|. On the other hand, the relation among the valence mass coefficients, Eq. (HOl) . 



is obeyed by the expressions for these coefficients found in Ref. 



33| . This occurs because 



the contributions of the terms proportional to and ps in Eqs. ( IHTl) and (l82l) . which were 
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not considered in Ref. [33j, are proportional to rrix + my and automatically obey Eq. (HOj) . 

Note, finally, that the relation in Eq. fHOl) is almost certain to be violated at next order 
in HQET. This is because the contributions from operators like the ki term in Eq. (ISTj) 
will affect the B and the B* differently at OII/itlq), destroying the cancellation that made 
Eq. (HO]) possible. 

V. FINITE VOLUME EFFECTS 

In a finite volume, we must replace the integrals in Eqs. fIBSp through fIBSp by discrete 
momentum sums. We assume that the time direction is large enough to be considered 
infinite (this is the case in MILC simulations), and that each of the spatial lengths has 
(dimensionful) size L. 



The correction to Eq. (IBSp is given explicitly in Ref. 12(]. In finite volume, we need only 
make the replacement 

/i(m) ^ /{^(m) = /i(m) + m'^6i{mL) . (83) 
Here 6i is a sum over modified Bessel functions 

where r is a 3- vector with integer components, and r = |r |. 

Arndt and Lin 42| have worked out the finite volume correction to Eq. (]B6|) . In our 
notation, the function l2{m, A) is replaced by its finite volume form, I^{m, A), 

him, A) f^{m, A) = him, A) + 5h{m, A, L) , (85) 

where the correction 5h{m, A, L) is given simply in terms of the function JpY{m, A, L) 
defined in Eq. (44) of Ref. '^-.^^ 

6h{m,A,L) = -(47r)2 A Jfv("^, A, L) 



We have added the L argument to Jpv for consistency with our notation 
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with ujq = ^J + 



The asymptotic form of JFY{m, A, L) for large mL is useful for practical applications, 
where typically mL > 3, and often mL > 4 [8|, [Oj. Arndt and Lin have found |42| : 

1 



Jfv("^, A,L) 
A 



8nrL 



-rmL 



A 



(87) 



e^' > [1 - Erf (2)] + 
1 



1 



rmL 
+ 



rmL ^ 
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1 



2; 



+ _e(^^)[l-Erf(;2) ] 



+ + + e(^^) h - Erf (z)l 



rmLy 



where 



z = 



rmL 



(88) 



(89) 



Computing higher orders in the l/{mL) expansion is possible if greater precision is needed. 

Since the functions Ii{m) and J2(m, A) arise from the integral 2^{m,A) in Eq. (]B7p . as 
well as from Eqs. fIBSP and flB6p . which serve to define them, it is necessary to check that 
the finite volume corrections coming from Eq. (]B7P are just those given by Eqs. (183!) and 
( l85|) above. This is easily seen to be true in the rest frame of the heavy quark, in which 
we are working. It is a consequence of the facts that: (1) in the rest frame, only the fi = 
component of X^{m,A) is non-zero, and (2) the integral over dq*^ is unaffected by finite 
volume, since we assume large time-extent of the lattices. The finite volume integral then 
splits into I^{m) and /|^(m. A) pieces, just as in infinite volume. 

Finally, we have to examine the finite volume corrections to the integral J^'^ , Eq. fIBSp . 
Since the function J2{m, A) does not enter our final results, we need only evaluate 



VuVfj.jq'^q" 



(27r)4 iy-q - A + ie){q^ - m^ + ie) 
d'^q —ic^ 
(27r)4 {vq - A + ie)(g2 - m? + ie) 

Ji(m,A) , 



(4vr) 



(90) 



where q is the spatial 3-vector part of q^. In the last line, the arrow refers to the fact 
that the function Ji arises after regularization and renormalization of the integral. A useful 
regulator in the present context is given by the insertion of a factor of exp(— cUg/Ag), where 
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Aq is a cutoff. After performing the contour integral over 



J 



(27r)3 2ujg{ujg + A) 

d^q 1 f d^q A f d^q A^ — 



(27r)3 2 7 (27r)3 2cug ' 7 (27r)3 2cjg(cj, + A) ' ''^^^ 
The first term is a pure divergence with no m or A dependence. It is thus the same in 



finite volume or infinite volume 12|. The correction to the middle term is proportional to 
the correction to Ji, since the same integral appears after performing the g° integration in 
Eq. ( ]B5I) . Similarly, the integral in the third term is proportional to that arising from the 
g° integration in Eq. ( ]B6I) . and the correction is therefore already known. We have 

Ji(m, A) ^ J["(m, A) = Ji(m, A) + 6Ji{m, A, L) , (92) 

where 

TTV^ — A^ Tn^ 
SMm, A, L) = 5h{m, A, L) - — 5i{mL) (93) 

The correction to Jf"*^, Eq. flSTj) . is 

6Jr\m,A,L)=6Mm,A,L) + ^^JMm,0,L) , (94) 

where Jpyim, 0, L) is the same as Eq. fl87|) with ^ = 1. 

With the expressions in this section, it is straightforward to incorporate the corrections 
to Ji, I2, and Ji numerically into fits to finite- volume lattice data. 

VI. CONCLUSIONS 

We have presented the NLO expressions in partially quenched SXPT for the form factors 
associated with B P^y semileptonic decays, for both infinite and finite volume. Using 
a quark flow analysis, we have obtained these results by generalizing the NLO PQXPT 



expressions calculated in the continuum in Ref. [331]. The main subtlety in applying this 
technique is due to the appearance of taste matrices inside the Feynman diagrams, since 
non-trivial signs can arise from the anticommutation relations of the taste generators. We 
have shown that these signs can be accounted for by a careful analysis of the relevant quark 
flow diagrams. 

The SXPT expressions are generally necessary for performing chiral fits to lattice simu- 
lations where staggered light quarks are used. For simpler quantities than the form factors, 
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SXPT has been seen to be essential [10|, l20| in order to get reliable extrapolations both to 
the continuum limit and to the physical quark mass values. For form factors, the lattice 
data in Ref. jo!] was not yet sufficiently precise for the SXPT expressions to be required (over 
continuum forms) for acceptable fits. However, we expect that the forms derived here will 
become more and more important as the lattice data improves. 

Our results are valid to lowest order in HQET; in general, we neglect I/tub corrections. 
We do however include the B*-B splitting A* on internal B* lines that are not in loops. 
This prescription allows the form factor fp to have the physical pole structure. Our 



treatment of the B*-B splitting is similar, but not identical, to that of Refs. [33|, |40| . Unlike 
those authors, we iterate self-energy contributions, namely Fig. M^a) and the effect of the 
one- loop mass shift of the B, to all orders. This seems to us to be a natural choice, and 
also makes the one-loop corrections better behaved. Indeed, with the values of light quark 
masses and momenta typically used in staggered simulations [s, Q], the one-loop B mass 
shift can dominate other one-loop corrections, so summing such self-energy contributions 
to all orders seems entirely appropriate. The final answers are then expressed in terms of 
the splitting A*^ = Mb* — Mb^ ■ In fitting lattice data, we suggest using the actual lattice 
values of this mass difference (at the simulated light quark mass values and lattice spacings), 
rather than applying a one-loop formula for the mass shifts. 

Our primary results for the staggered, partially quenched case with three non-degenerate 
sea quarks are found in Sec. IIV A[ The form factor (also known as /y) is given by 
Eq. ( IMl) in terms of quantities defined in Eqs. ( 135|) . ( 152|) and ( l53i) . as well as the wave 
function renormalization factors SZp^^ and SZb., that are listed in Eqs. (IClj) and (1C2I1 of 
Appendix O Similarly, the form factor fp (also known as f±), is given by Eq. (136|) in terms 
of quantities defined in Eqs. fl3T|) . (|37|) . fl38|) . and fl5^ through fl57j) . as well as the wave 
function renormalization factors. We have also found a single relation, Eq. (140!) . among the 
parameters that control the analytic valence mass dependence. While this relation is also 
satisfied by the parameters written down in Ref. [33|, it is important to know that it persists 
even in the presence of the complete NLO forms of the Lagrangian and current. 

Appropriate limits of our expressions can be taken for various relevant cases, including 
the case of full (unquenched) staggered QCD (in Sec. IIVB|) and the case of continuum 
PQXPT with non-degenerate sea quark masses [Eqs. fl63|) . fl64l) . flC3l) . and flC4|) ]. Despite 
the fact that the latter are continuum results, they have not, to our knowledge, appeared 
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in the literature before. Finally, our expressions can be corrected for finite volume effects 
using the results of Sec. IVl 
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APPENDIX A: FEYNMAN RULES 

In this appendix we list the SXPT propagators and (some of) the vertices in Minkowski 



space 



19| . as well as the corresponding continuum versions. 



In SXPT, the propagators for the heavy-light mesons are 

Here a, h indicate the flavor-taste of the light quarks, and A* is the B*-B splitting in the 
chiral limit, which we often neglect since we work to leading order in HQET. 
The 33*71 vertex is: 

J {BlB;,-3;i3,) d>^<^,,, (A3) 

where repeated indices are summed. Other needed vertices come from the expansion of the 
LO current, Eq. (|T8l) . We have: 



J^'^ = k3*/ (^5ac - ^$ab$.c + ■■■)+ KV^3a (^^<I>., + ■ ■ ■ ^ , (A4) 

where repeated indices are again summed and ■ ■ ■ represents terms involving higher numbers 
of pions, as well as contributions from the axial vector part of the current, which are not 
relevant to the form factors. 

If desired, each fiavor-taste index can be replaced by a pair of indices representing flavor 
and taste separately. We use Latin indices in the middle of the alphabet (i, j, . . . ) as pure 
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flavor indices, which take on the values 1, 2, . . . , A'sea in full QCD. Greek indices at the 
beginning of the alphabet 7, . . . ) are used for quark taste indices, running from 1 to 

4. Thus we can replace a ^ ia and write, for example. 



As in Refs. 



13|, 



2(vk + ie) 



(A5) 



19| . pion propagators are treated most easily by dividing them into 



connected and disconnected pieces, where the disconnected parts come from insertion (and 
iteration) of the hairpin vertices. The connected propagators are 



k J conn 



- mf- .^ + ie 



where S is one of the 16 meson tastes [as defined after Eq. ([5]) 
mass of a taste-H meson composed of quarks of flavor i and j: 

mfj^ = jiirrii + rrij) + a^Ac. 



(A6) 

and mij^s is the tree-level 



(A7) 



Here A= is the taste splitting, which can be expressed in terms of Ci, C3, C4 and Cq in 
Eq. (JT6II jisl . There is a residual 5*0(4) taste symmetry H] at this order, implying that 
the mesons within a given taste multiplet (P, V , T, A, or J) are degenerate in mass. We 
therefore usually use the multiplet label to represent the splittings. 

Since the heavy-light propagators are most simply written with flavor-taste indices, as in 
Eqs. ( IA1[) and (1A2[) . it is convenient to rewrite Eq. (]A6[) in flavor-taste notation also: 

iSii'6jj'T^i^Ti^,^, 



^ J conn 1^ J conn p — iil^^ ~r ^6 



(A8) 



where are the 16 taste generators, Eq. IQ. 

For flavor- charged pions {i 7^ j), the complete propagators are just the connected propa- 
gators in Eq. flA6P or flASp . However, for flavor-neutral pions {i = j), there are disconnected 
contributions coming from one or more hairpin insertions. At LO, these appear only for taste 
singlet, vector, or axial-vector pions. Denoting the Minkowski hairpin vertices as —iS'^, we 



have 



13|: 



6L 



T-E € {^fi} (taste vector); 
Ts G {^/is} (taste axial- vector) ; 
Aml/S, Th = ^/ (taste singlet); 
0, 7= £ {^fiuj^b} (taste tensor or pseudoscalar) 



(A9) 
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with 



16 

^V{A) = -piC2V{A) - C5V(A)) ■ (AlO) 



The disconnected pion propagator is then 

{^^^^%'}^JP^ = ' (All) 



where 

For concreteness we have assumed that there are three sea-quark flavors: m, c?, and s; the 
generahzation to A^sea flavors is immediate. Here m{/= = muu,s is the mass of a taste-S 
pion made from a -u and a u quark, neglecting hairpin mixing (and similarly for m£),H and 
ms^^), "ni^o H, m^ H; and m^/ = are the mass eigenvalues after mixing is included, and the ie 
terms have been left implicit. When specifying the particular member of a taste multiplet 
appearing in the disconnected propagator is unnecessary, we abuse this notation slightly 
following Eq. (1A7P and refer to ^^J^j/j/, 1^ flavor-taste notation we have: 

|$a6$6'a4 = {$^a,i/3$i'/3^^'a'| (p) = ^ij^jV V TJjTI^,©! ,,,, (Al3) 

k J disc I. J disc ' ' 



For comparison, we now describe the continuum versions of the Feynman rules [18j. Since 
taste violations do not appear in Cul, Eq. (fT7|) . the continuum-theory version of Eqs. (lAip 
and flA2l) are unchanged except that flavor-taste indices are replaced by pure flavor indices 

{b;.B:J}(*) = "^l^'r;,^)""' l»ntin-m]. (A15) 

Similarly, the continuum BB*tx [isj] and current vertices are identical to those in SXPT, 
aside from the redeflnition of the indices and a factor of 2 for each $^5 fleld due to the 
non-standard normalization of the generators in the SXPT case, Eq. ([6]). The continuum 
version of Eq. (1A3[) is 

2^ (5;! B, - Bj B;^i d^^y. [continuum]; (A16) 
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while the continuum version of Eq. (lA4p is 

jCo' = f^B*/ (^Sik - ^^ii^ik + ■ ■ ■ ^ + fivf'Be + ■ ■ • ^ [continuum]. (A17) 

Because of taste-violations in the SXPT pion sector, the differences between the propaga- 
tors Eqs. flA6p . flAllI) and flA12p and their continuum versions are slightly less trivial. The 



U.,<l>j'A ip) = o ' V'' [continuum], (A18) 

I- J conn U — nl-- + it 



continuum connected propagator is 

■p^ - mfj + le 
with 

rn^j = fx(mi + rrij) [continuum]. (A19) 
The continuum disconnected propagator is 



<^$ij<l>j'i/i (p) = Sij6j'i'Vii^i>i> [continuum], (A20) 

L J disc 



where [13 1 



V,i,i.e = -^S'-^ 2T71 TTTi fvTl fr [continuum], (A21) 

with now 6' = itlq/S. 

Note the difference in normalization between S' and the SXPT taste-singlet hairpin, 5j, 
Eq. flA9p . This arises from the fact that itlq/S is defined to be the strength of the hairpin 
vertex when one has a single species of quark on each side of the vertex [l^ . In the staggered 
case, each normalized taste-singlet field is made out of four species (tastes), for example 
(p^ = I (011 + 022 + 033 + 044), whcrc is flavor neutral, and only taste indices are shown. In 
the disconnected propagator of two such fields, there are 16 terms, and a factor of (1/2)^ from 
the normalization, so there is an overall factor of 4 relative to a single-species disconnected 
propagator, such as that of 0ii with 022- At one loop, the "external" fields in this propagator 
are always valence fields, so the normalization issue has nothing directly to do with the 
fourth root trick for staggered sea quarks. (The normalization is in fact compensated by the 
extra factors of 2 in the continuum vertices.) The rooting does however affect the rjj mass 
that appears in denominator of Eq. (1A21I) . which comes from iterations of the hairpin and 
therefore involves sea quarks. The end result is that m^, j ^ Nseafnl/S (for large mo), rather 
than pa 4A^sea^o/35 the value in the unrooted theory [13|. In the continuum, we also have 
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APPENDIX B: INTEGRALS 



Here we co 
factors 



19 



33| 



lect the integrals needed in evaluating the diagrams for the semileptonic form 



The disconnected propagators can be written as a sum of single or double poles using 
the (Euclidean) residue functions introduced in Ref. or their Minkowski-space versions. 
We define {m} = {mi, m2, ■ ■ ■ , rrin} as the set of masses that appear in the denominator of 
Eq. (]A12p . and {^} = {/ii, /X2, • • • , /Wfc} as the numerator set of masses. Then, for n > A; and 
all masses distinct, we have: 



n,=i(g^-m2 + ie) 



where the Minkowski space residues Rf''^'^ are given by 



— m, + le 



(Bl) 



nliK-/^?) 



If there is one double pole term for = m'j (where mi G {tti}), then 



(B2) 



(g2 ~mj + ie) n"=i(?^ - ^"j + 

d 

dm?n 



(B3) 



In the end we want to write the results in terms of the Euclidean-space residues R^p''\ 
because they are ones we have used previously [isl, Q- In Euclidean space the sign of each 
factor in Eq. ( ]B2[) is changed. We therefore have 

nli(/^? 



n 



I2 



-^3 



The integrals needed for the form factors are 
f d'^q i 1 



17 



33|) 



/i 



4~d 



d'^q i 



him) 



{27rY (v-q - A + ie)(g2 - m'^ + ie) (Att)^ A 



him, A) 



A-d 



iq^ 



J- 



111/ 



(27r)'^ iv-q — A + ie)iq^ — m^ + ie) 
d'^q 



iq^q" 



(47r)2 
A 



[him, A) + him)] 



(B4) 

(B5) 
(B6) 
(B7) 



(27r)'^ (t;-g- A + ie)(g2 -m^ + ie) (47r) 



[Ji(m, A)g^"' + J2(m, A)v^v''] 



(B8) 
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where the arrows represent the fact that the r.h.s. of these expressions have already been 
renormahzed (unhke the corresponding equations in Ref. [33| ) . 



APPENDIX C: WAVEFUNCTION RENORMALIZATION FACTORS 



The one loop chiral corrections to the wave function renormalization factors Zb and Zp 
are are [13I, Q 



bZi 



3(47r 
1 

+ :t 



{ 16 |J 



+ 2 i?;^'^l(A^r-);/.f)/,K,,) 



I 2 c-/ 



dm 



V, 



yy 



2 ^r'[^f^''-^^f)^M..) 



+ 



(CI) 



bZ 



(4vr/)^ 



- E 



9 



3 dm\ T 



(3,:r). , (3) 



d 



dm\ y L 



RfHMr':i^f)Hm,y) +[V 



where / runs over the sea quarks (m, rf, s). 



A] , (C2) 
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For the continuum result in partially quenched XPT, we can simply set a = and ignore 
taste splittings. In the 1 + 1 + 1 case, we get 



1 

+ 3 



.jeX(3,x) ^ 

d 



+ 2 i??'''(-^^''''^;/^^'ViK) 



(C3) 



-Da: 



(C4) 



Returning to a 7^ 0, and taking the valence quark masses to be m-^ = rUy — = m^, we 
have the 2 + 1 full QCD pion result in SXPT: 



-Ii{m^y) + 



+ (-4a25' 



(m2^ -m2^)(m2^ -^^^)' 



+ 



(C5) 



Taking the valence quark masses to be rrix = m„ = rrid and niy = nis gives the 2 + 1 full 
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QCD kaon result: 



3(47r 



1 3 

--/i(m^J + -/i(m^J - /i(m5j 



fm?. — rni ) (171% — rni ) iim?, — rnl ) 
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Setting m^ = m^ = in Eq. flC2p results in the 2 + 1 full QCD result for the B 
wavefunction renormalization: 
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Finally, putting m^ = mg and m„ = md in Eq. ( ]C2[) . we obtain the full QCD Bs renormal- 
ization factor in the 2 + 1 case: 
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TABLE I: Connecting the one-loop diagrams from Ref. [3i 
column) . 



(left column) and this paper (right 



xiei. [o^j 
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ig- l^aj 
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(13) 
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(14) 


Fig.Hb) 




FIG. 1: Example of a connected one- loop form factor diagram at (a) the meson level and (b) the 
quark level. For the meson diagram, the double line is a heavy-light meson while the single line is 
a pion. For the quark-level diagram, the solid line is a heavy quark and the dashed line is a light 
quark. The internal sea quark loop is required by the (quark-flow) connected pion propagator; 
purely valence diagrams are only possible with a disconnected pion propagator. Therefore this 
diagram gives rise to a factor of A'sea in the degenerate case. 
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(«) 



(6) 



FIG. 2: Example of a disconnected one-loop form factor diagram at (a) the meson level and (b) 
the quark level. The cross in the meson diagram represents the two-point interactions in XPT, 
and is represented by the "hairpin" in the quark-level diagram. There are no factors of A'^sea but 
instead factors of 1/A^sea coming from the decoupling of the ry'. 



(a) 



ib) 



FIG. 3: Tree level diagrams for (a) and (b) fp. The double line is the heavy-light meson and 
the single line is the pion. 
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(a) 



(&) 



FIG. 4: One-loop diagrams. The internal light meson lines may in general be connected or 
disconnected: possible hairpin insertions are not shown explicitly. 




(a) 







1 1 




1 





FIG. 5: One-loop fp diagrams. The internal light meson lines may in general be connected or 
disconnected: possible hairpin insertions are not shown explicitly. 
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FIG. 6: The quark- flow diagram for Fig. [5]^b), omitting the heavy quark line for clarity. The 
mesons in the loop are X and Y mesons, flavor-neutral mesons made up of x and y quarks. Note 
that even though only a single hairpin insertion is shown explicitly, the figure should be interpreted 
as representing all diagrams with one or more hairpins. 
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FIG. 7: Possible quark- flow diagrams for Fig. [5^c) with a disconnected meson propagator in the 
loop. The solid rectangle encloses the 5-point vertex of Fig. [5]^c). The heavy quark line has been 
omitted for clarity. A "reflected" version of diagram (b), with the outgoing pion on the other side 
of the vertex, is also possible. 
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(a) (6) 

FIG. 8: Possible quark-flow diagrams for Fig. [SI^c) with a connected meson propagator in the loop. 
The solid rectangle encloses the 5-point vertex of Fig. [5l|c). The heavy quark line has been omitted 
for clarity. Since we have assumed that x and y are different flavors, diagram (a) cannot occur in 
our case. Diagram (b) can occur, as can a "reflected" version with the outgoing pion on the other 
side of the vertex. 
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